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Abstract 

Single channel measurements play a minor role in today physics, but they are some- 
times unavoidable. Comparing to multichannel measurements, there is distribution 
of measurement time to be chosen in an experiment design. A method to optimize 
distribution of measurement time is given, where optimal distribution minimizes 
standard deviation of a selected fit parameter. As an example, the method is ap- 
plied to electron spectroscopy experiments. 
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1 Introduction 



Single channel measurements offer a unique challenge. Regarding a setup of 
such experiment, there is one more characteristics to be fixed in advance during 
the experiment setup, compared to multichannel measurements. This is the 
distribution of measurement time. 



Let there be a single channel experiment in its design phase. A spectrum to be 
recorded by the experiment is dependent on several parameters. Experimenters 
are often interested in one of them. Let the mechanical design of experiment 
be fixed, and let the measurement points be chosen. So, it is known how 
to measure, and where to measure. However, it is unknown for how long to 
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measure in the particular measurement points, which we call a distribution of 
measurement time. 

Then the task is to distribute the measurement time into the measurement 
points in the way, that minimizes the standard deviation of the parameter, 
in which the experimenters are interested. Up to our knowledge, there are 
no studies available concerning this aspect. The point is, we can not set the 
measurement times for all the points simultaneously (it is numerically not 
feasible to minimize a function of tens or hundreds of parameters). Instead, 
we can fix the measurement time in points one by one, approaching the optimal 
distribution of measurement time in an iterative way. 



2 Method 



Let {Ei}, i = l...n be a fixed set of measurement points, where n is their 
number and r to t is the total measurement time. Further, let {T{\ be the initial 
time distribution chosen intuitively, e.g. uniform distribution, meaning that 
equal time is spent in each point. Let the initial time distribution result in 
an initial standard deviation of the selected fit parameter cr init = a({Ti}) that 
was obtained by simulations (creation of pseudo-experimental spectra and 
their evaluation). 

Now, we minimize standard deviation of the selected parameter varying the 
time Ti in the first measurement point only. (This means that we simulate 
spectra corresponding to various 7\ and evaluate them to find the value of Ti 
supplying the minimum standard deviation of the selected fit parameter a. 
The measurement times in all the rest points are scaled by the same factor to 
keep the total measurement time equal to r tot . As the result we obtain the time 
distribution {T^} and the corresponding standard deviation of the selected 
fit parameter = cr({T/^}). So, the measurement time in the first point is 
T±~\ and the measurement time in each of the rest points is equal to 

Ttot - 1\ 



i = 2 ... n. Then we do the same for all the other measurement points (starting 
always with the {Tj} distribution). So, we get n time distributions {T- k ^} , k = 
1 ...n and corresponding sigmas a^ k \ (For the sake of clarity, we repeat that 
o~( fe ) are standard deviations of the same parameter but obtained from spectra 
corresponding to different time distributions T- h \) 

Now we construct the best final time distribution from the set {T- k ^} of the 
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partial ones. There are several possibilities how to do it. Experimentally, it 
turned out that a good way is a weighted sum 

1 n 

^ = -E^f } - (2) 

n k=i 



The weight factors uik may be estimated as: 



where s is chosen to minimize a({T(}), i.e., we simulate spectra corresponding 
to various s and evaluate them to find the value of s minimizing the standard 
deviation a. 

Having the new time distribution {T(} we renormalize it with respect to the 
Ttot- (Note, that the weight factors in eq. ([3]) do not keep the J2T( = r to t-) 
Then, we use it instead to the original {Tj} to repeat the whole process in an 
iterative way. 

Based on our extensive simulations, this seems to be the best way how to com- 
bine the partial time distributions {2^ } together, guaranteeing convergence, 
numerical stability, and reasonable speed. The partial time distributions that 
significantly improved the standard deviation of the selected parameter are 
favored, and the convergence speed is enhanced by the choice of the s factor. 
(Note, that eq. ([3]) is well defined because a in i t > a^ k ' for any k.) 

We also tried several other ways how to merge the partial time distributions 
{T^} together. In particular, following choice of the final time distribution: 

T; = Tf * I™ (4) 



looked feasible and promised fast convergence to the optimal time distribution. 
Nevertheless, it proved to be highly numerically unstable, and it oscillated 
practically always. Finally, if the numerical stability is favored, we tested that 
setting all the u k in eq. ([2]) equally to one was a good choice. However, the 
convergence is about 5 times slower. Practically, we implement the method 
using a simplex minimizer provided by Minuit2 library wrapped inside the 
ROOT framework [TJ. 
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Fig. 1. An example of the optimal distribution of measurement time. The time dis- 
tribution minimizes the standard deviation of Gaussian line position fitting four 
parameters: amplitude, background, position, and width. The initial parameter val- 
ues are given in text. For convenience the Gaussian line is shown as well. 



3 Examples 



To test the method we have chosen an idealistic Gaussian line on a constant 
background as measured by a differential spectrometer: 



\2 N 



G(E\A, B, E Q , a) = A • exp ( 1 + B . (5) 



The line is described by four fit parameters: amplitude A, background B, 
position E and width a. The initial parameter values were chosen as follows: 
the amplitude of 10 Hz, the background of 10 Hz, the line position of 25 eV, 
and the width (sigma) of 5 eV. The distribution of measurement points was 
chosen to be an uniform one from eV up to 50 eV with a step of 1 eV. 

We studied two cases; optimal time distribution with respect to the standard 
deviation of the Gaussian line position and that with respect to the line width. 
As for the line position, the resulting optimum time distribution is given in 
fig. [TJ It turned out, surprisingly, that the most of the relevant information 
for this case is concentrated in four points only. The standard deviation of 
the line position improved by a factor of 2.0 with respect to the uniform time 
distribution. When optimizing the line width standard deviation, we reached 
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Fig. 2. An example of the optimal distribution of measurement time. The time 
distribution minimizes the standard deviation of Gaussian line width fitting four 
parameters: amplitude, background, position, and width. The initial parameter val- 
ues are given in text. For convenience the Gaussian line is shown as well. 

the time distribution depicted in fig. [2j Here, the particular standard deviation 
improved by a factor of 1.4. (In both cases, the optimization procedure needed 
100 iterations.) 

As the second example, we have chosen an electron spectroscopy experiment 
to determine the neutrino mass, in particular, the KATRIN experiment [2113] 
aiming at neutrino mass sensitivity of 0.2eV/c 2 . The necessity of the single 
channel measurement originates from the intrinsic property of the applied 
spectrometer type, that is the only one exhibiting simultaneously required 
energy resolution and luminosity. There are four fit parameters in the experi- 
ment: the neutrino mass, beta spectrum endpoint, amplitude and background. 
Integrated beta spectrum as proposed to be measured by the KATRIN exper- 
iment is given in [3], as well as the standard distribution of the measurement 
points, which is a non-equidistant one. 

First, let the neutrino mass be the selected parameter. Then, starting with 
the uniform time distribution, and optimizing the time distribution in 16 iter- 
ations, the simulation resulted in the optimal time distribution shown in fig. 
|3j The standard deviation of the neutrino mass was improved by factor 1.18 
compared to the uniform time distribution. The same improvement could be 
achieved by a prolongation of the total measurement time by factor of 1.9. 
(Standard deviation of the neutrino mass scales with fourth root of the total 
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Fig. 3. An example of the optimal distribution of measurement time, which could be 
used by the KATRIN experiment. The time distribution minimizes the standard de- 
viation of the neutrino mass when fitting four parameters: amplitude, background, 
the tritium /3-spectrum endpoint, and the neutrino mass. Note, the initial free neu- 
trino mass of eV, and the initial free endpoint of 18 575 eV was chosen to run the 
simulations. 

measurement time, since the neutrino mass squared enters the form of tritium 
beta spectrum.) 



4 Other applications 

The method can also be effectively used to show what energy region of the 
measured spectrum affects mostly a standard deviation of a selected parame- 
ter. As an example, we focused on the endpoint of beta spectrum. Since the 
beta spectrum endpoint and the neutrino mass are strongly correlated [2|3] . 
we fixed the neutrino mass to OeV/c 2 . Then, the optimal time distribution 
minimizing the standard deviation of the endpoint was derived. It is shown in 

fig.H 

The optimal distribution of measurement time was found following the method 
described in the above paragraphs using the same standard distribution of the 
measurement points as in the previous example. Here, it was the standard 
deviation of the tritium endpoint we were focused on, not the standard devi- 
ation of the neutrino mass. Again, starting with a uniform time distribution 
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Fig. 4. An optimal time distribution with respect to the minimal standard deviation 
of the endpoint energy assuming a fixed and known neutrino mass, and fitting am- 
plitude, background, and the endpoint energy. Compare with the time distribution 
in fig. [3] optimized with respect to the standard deviation of the neutrino mass. 

10 iterations were performed. The method sets longer measurement times in 
the measured points that are sensitive to the beta spectrum endpoint. 

We would like to note that the same method can be used to suppress sys- 
tematics, i.e., to find a region of the spectrum, which is the most sensitive to 
systematics, and then to exclude the region from the set of the measurement 
points. 



5 Discussion 

A method to distribute the measurement time in single channel measurements 
into the measurement points to minimize the standard deviation of the selected 
fit parameter was offered and demonstrated on two examples. It worked the 
desired way. Even more, the method proved useful to show what energy region 
of a measured spectrum is sensitive to a selected parameter. 

In place of summary, we would like to emphasize, that the method is mathe- 
matical. And so are its results. Although the method assumed a fixed distri- 
bution of measurement points, this assumption should be reconsidered if nec- 
essary. E.g., if the optimal distribution of measurement time aggregates most 
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of the measurement time into few measurement points, then a more dense dis- 
tribution of points in the particular measurement region is appropriate. This 
actually happened in the KATRIN example. The case of measurement points 
with negligible measurement time can be treated in a similar way. An extraor- 
dinary care should be paid to when regions sensitive to the selected 
parameter overlap regions sensitive to systematics. 
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